to appear in The Astrophysical Journal 



STIS Echelle Observations of the Seyfert Galaxy NGC 4151: 
Physical Conditions in the Ultraviolet Absorbers^ 

S. B. Kraemer^, D. M. Crenshaw^ J. B. Hutchings^ I.M. George^'^ A.C. Danks^, T.R. 
^ _ GuU^ M.E. Kaiser^ C.H. Nelson^ & D. Weistrop^ 

O 
(N 

55 ■ ABSTRACT 



^ . We have examined the physical conditions in intrinsic UV-absorbing gas in 

the Seyfert galaxy NGC 4151, using echelle spectra obtained with the Space 
Telescope Imaging Spectrograph (STIS) on the Hubble Space Telescope (HST) 
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. on 1999 July 19. We confirm the presence of the kinematic components detected 

^ in earlier Goddard High Resolution Spectrograph (GHRS) observations, all of 

which appear to be outflowing from the nucleus, as well as a new broad absorption 
Q . feature at a radial velocity of -1680 km s~ . The UV continuum of NGC 4151 

^ ! was a factor of about 4 lower than in observations taken over the previous two 
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years, and we argue the changes in the column density of the low ionization 
absorption lines associated with the broad component at -490 km s~^ reflect the 
decrease in the ionizing flux. Most of the strong absorption lines (e.g., N V, C IV, 
Si IV, etc.) from this component are saturated, but show substantial residual 
flux in their cores, indicating that the absorber does not fully cover the source of 
emission. Our interpretation is that the unocculted light is due to scattering by 
free electrons from an extended region, which reflects continuum, emission lines, 
and absorption lines. For the first time in such a study, we have been able to 
constrain the densities for this kinematic component and several others based on 
the strength of absorption lines from metastable states of C III and Fe II, and/or 
the ratios of ground and fine structure lines of O I, C II, and Si II. We have 
generated a set of photoionization models which not only successfully match the 
ionic column densities for each component during the present low flux state, but 
also those seen in previous high flux states with the GHRS and STIS, conflrming 
that the absorbers are photoionized and respond to the changes in the continuum 
flux. Based on the model parameters (ionization parameter and density), we 
have been able to map the relative radial positions of the absorbers. We find 
that the absorbing gas decreases in density with distance. Finally, none of the 
UV absorbers is of sufficiently large column density or high enough ionization 
state to account for the observed X-ray absorption, while the scatterer is too 
highly ionized. Hence, the X-ray absorption must arise in a separate component 
of circumnuclear gas. 

Subject headings: galaxies: individual (NGC 4151) - galaxies: Seyfert - ultravi- 
olet: galaxies 

1. Introduction 

NGC 4151 (cz = 995 km s^^; determined from 21 cm observations, de Vaucouleurs et al. 
[1991]) is the first Seyfert galaxy known to show intrinsic absorption that could be attributed 
to the active nucleus. Ultraviolet observations of NGC 4151 by the lUE (Boksenberg ct 
al. 1978) and subsequent far-ultraviolet observations by the Hopkins Ultraviolet Telescope 
{HUT, Kriss et al. 1992) revealed a number of absorption lines from species that span a wide 
range in ionization potential (e.g., O I to O VI), as well as fine-structure and metastable 
absorption lines. The intrinsic UV absorption was found to be variable in ionic column 
density, but no variations in radial velocities were detected (Bromage et al. 1985). The 
ionization state of the UV absorbers are generally well correlated with the UV continuum 
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flux level; the equivalent widths of low ionization species are negatively correlated with the 
flux level, while the those of the high ionization species are positively correlated (Bromage 
et al. 1985). More recent observations of NGC 4151 were obtained with the Goddard High 
Resolution Spectrograph (GHRS) at high spectral resolution (~15 km s~^) over limited 
wavelength regions by Weymann et al. (1997). The GHRS spectra revealed that the C IV 
and Mg II absorption lines, detected in six major kinematic components, were remarkably 
stable over the time period 1992 - 1996, with the notable exception of a transient feature 
which appeared on the blue wing of the broad C IV emission. However, this feature, which 
was discussed in our previous paper (Crenshaw et al. 2000b; hereafter Paper I), turns out 
to be due to Si II, rather than a high velocity C IV component, as suggested by Weymann 
et al. (1997), 

X-ray spectra of NGC 4151 reveal the presence of a large intrinsic column of atomic gas 
(see, e.g., Barr et al. 1977; Holt et al. 1980). Based on ASCA spectra (see George et al. 
1998, and references therein), the 0.6 - 10 keV continuum of NGC 4151 can be fit as a fiat 
power law (a = 0.3 - 0.7), with an absorber of column density ~ 1022-23 cm~^, which is 
most likely ionized, as suggested by Yaqoob, Warwick, & Pounds (1989). NGC 4151 is highly 
variable in X-rays (cf. Perola et al. 1982, 1986), and the X-ray absorber appears to respond 
to changes in the 2-10 keV flux, although the changes in total column and ionization state 
are not simply correlated with the X-ray luminosity (George et al. 1998). For example. 
Weaver et al. (1994) found evidence that the ionization state and column density of the 
absorber were both higher when the source was fainter, which may be more easily explained 
by transverse motion. 

In the Hubble Space Telescope {HST) survey of intrinsic UV absorption lines in Seyfert 
galaxies by Crenshaw et al. (1999), NGC 4151 stood out as an unusual object. Of the 
ten Seyfert galaxies in this sample with intrinsic absorption, all showed C IV and N V 
absorption, but only NGC 4151 showed Mg II absorption (although, recently, the presence 
of Mg II absorption has been confirmed in the Seyfert 1 galaxy NGC 3227; Kraemer et 
al. 2000b). NGC 4151 is also unusual in the sense that it is the only active galaxy known 
to show metastable C HI* A1175 absorption, which is indicative of relatively high electron 
densities (~ 10^ cm~^, Bromage et al. 1985). 

We obtained STIS echelle spectra of the UV spectrum to study the intrinsic absorption 
in detail. These observations are part of a long-term project on NGC 4151 by members of the 
Instrument Definition Team (IDT) of the Space Telescope Imaging Spectrograph (STIS) on 
HST. In Paper I, we concentrated on the variability of the absorption features in the spectral 
region surrounding the broad C IV emission. In this paper, we present results on the entire 
UV spectrum and the kinematics and physical conditions in the intrinsic absorbers. 
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2. Observational Results 

We observed the nucleus of NGC 4151 with a 0".2 x 0".2 aperture and the STIS E140M 
and E230M gratings to obtain a spectral coverage of 1150 - 3100 A at a velocity resolution of 
7-10 km s^^ (the ''STIS2" dataset). In Paper I, we give the details of the STIS observations 
and the data reduction procedures. We obtained the STIS echelle spectra on 1999 July 19, 
when the continuum and broad emission lines were in a low state compared to previous 
STIS and GHRS observations. We found that the UV continuum decreased by a factor 
of ~4 over the previous two years, which resulted in a dramatic increase in the column 
densities of the broad components of the low-ionization absorption lines (e.g.. Si II, Fe II, 
and Al II). In addition to the absorption lines seen in lUE spectra during previous low 
states (Bromage et al. 1985), we identified numerous broad Fe II absorption lines that arise 
from metastable levels as high as 4.6 eV above the ground state. We concluded that the 
low-ionization absorption lines are very sensitive to changes in the ionizing continuum, and 
that the extreme variability in these lines could explain a couple of puzzling aspects of the 
UV spectrum of NGC 4151: the "transient absorption feature" in the blue wing of the G IV 
emission feature (Weymann et al. 1997) and the "sateUite emission fines" (Ulrich et al. 1985; 
Glavel et al. 1987). 

2.1. Line Identification 

Figure 1 (a - f) shows the STIS echelle spectra over the entire UV bandpass. The 
broad absorption features (FWHM = 435 km s~^) appear at a hefiocentric redshift of 0.0017 
(velocity centroid = —490 km s~^ with respect to the nucleus), and correspond to a blend of 
kinematic components D and E in Weymann et al. (1997) We give our line identifications 
of the broad absorption features above the spectrum shown in Figure 1. Most of the line 
identifications at the top of Figure 1 agree with those claimed in lUE high-dispersion spectra 
by Bromage et al. (1985), including numerous resonance lines, excited fine-structure lines 
of G II, O I, and Si II, and the metastable G III* A1175 blend. In addition we have firm 
identifications of Ni II lines with oscillator strengths / > 0.04 (Morton, York, & Jenkins 
1988). Broad Gr II and Mn II absorption may also be present, but these identifications 
are less certain and may be affected by the numerous Fe II absorption lines in the spectra 
(see below). We find no evidence in the STIS spectra for the metastable He II A1640 and 
metastable He I A2764, A2829, and A2945 fines claimed by Bromage et al. (1985), although 
the Fe II lines may be masking these features as well. 

Below the STIS echelle spectra in Figure 1, we plot all of the Fe II UV fines in the fists 
of Silvis & Bruhweiler (2000), which include multiplets from levels as high as 4.6 eV above 
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the ground state. Few regions in the spectra are unaffected by Fe II absorption; this can be 
seen by noting that local regions of the spectra with no Fe II multiplets plotted below them 
(and no other absorption features above) are much higher than surrounding regions (e.g., 
1341 - 1346 A, 1593 - 1602 A, 2648 - 2662 A, 3014 - 3044 A). It is also clear that multiplets 
from levels high above the ground state must contribute; for example, Fe II UV multiplet 
HU, arising from a level 4.1 eV above ground, is clearly present in the observed spectrum at 
1310 - 1314 A (although the shorter wavelength hues are blended with Si II* A1309.3) 

The narrow kinematic components of absorption identified by Weymann et al. (1997) 
are also seen in a number of ionic species. Figure 2 demonstrates this point with plots 
of the Si IV A1393.8 and S II A1253.8 absorption as a function of radial velocity (with 
respect to the systemic redshift of 0.00332). For the relatively high ionization fine of Si IV, 
components A, B (Galactic), C, D+E (the broad feature), F, and F' are evident (F and F' 
likely arise in the interstellar medium or halo of NGC 4151, Weymann et al. 1997). For the 
low ionization line of S II, components B, D+E, E', and F are present (E' is not seen in the 
high ionization lines). The presence or absence of kinematic components in different ions 
provide valuable information on the physical conditions in each component, which will be 
explored in this paper. An expanded atlas of the STIS echelle spectra with our identification 
of each component of absorption is given in Danks et al. (2000) . 



2.2. Saturation and Scattered Light 

There are a number of Galactic (component B) lines that are heavily saturated (e.g.. 
Si III A1206.5, Si II A1260.4, C II A1334.5, Mg II AA2796.3, 2803.5), as evidenced by their 
extremely small residual intensities. The troughs of these Galactic features all lie within 
1% of the nearby continuum fiux levels, indicating not only heavy saturation, but accurate 
subtraction of the instrumental scattered light as well (see Paper I) . Although the residual 
intensities in the troughs of the broad (D+E) absorption fines are much higher (>7% of 
the continuum in all cases), several clues indicate that most of the broad absorption lines 
are actually highly saturated, and the large residual intensities are due to unabsorbed (e.g., 
scattered) fight in NGG 4151. 

Galculation of the covering factor in the line of sight using the doublet method (Hamann 
et al. 1997; Grenshaw et al. 1999) yields values significantly less than unity in the troughs: 
0.91 ± 0.02 for N V AA1238.8, 1242.8 and G IVAA1548.2, 1550.8; 0.85 ± 0.02 for Si IV 
AA1393.8, 1402.8; and 0.62 ± 0.04 for Mg II AA2796.3, 2803.5. This indicates unabsorbed 
emission in the troughs that is not constant as a function of wavelength. In addition, the 
presence of a number of Fe II lines with small 5/ values (Silvas & Bruhweiler 2000), indicates 
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that Fe II lines with much higher values must be saturated. For example, in the 2900 - 
3100 A region in Figure If, the Fe II lines in multiplet 78 with moderate g/ values are nearly 
the same strength as the line with the largest value (at 2990 A), and even the line with the 
lowest gf value in this multiplet is relatively strong. Lines in other multiplets with much 
higher gf values are clearly saturated, although the residual intensities in their troughs can 
be quite large. Finally, the strengths of broad absorption lines with low oscillator strengths 
and/or low abundances (e.g., S II, Ni II) indicate that many of the other low-ionization 
absorption lines should be saturated. Thus we have strong evidence for heavy saturation of 
the D+E component in most of the absorption lines and a significant amount of UV light in 
the nucleus that is not occulted by these broad absorbers. 

There are indications that the spectrum of the unocculted light is not simple (e.g., a 
power-law continuum). A power- law continuum drawn through the majority of absorption 
troughs in Figure 1 would he well above the La, N V, Si IV, and C IV absorption troughs, 
which suggests a lower contribution in these regions (e.g., due to absorption features in 
the unocculted spectrum). On the other hand, this continuum would lie well below the 
level required to account for our measured low covering factor for the Mg II absorption. In 
addition, at the densities required for the presence of the Fe II mctastable absorption (> 10^ 
cm~^; Wampler, Chugai, & Pettijean 1995) or the C III metastablc absorption (> 10^ cm~^; 
Bromage et al. 1985), the Si II fine-structure line at 1533.4 A (Figure lb) should be twice 
as strong as the Si II resonance line at 1526.7 A (Flannery, Rybicki, & Sarazin 1980), and 
it does not appear to be so. The only simple solution to this dilemma is that both of these 
lines are saturated and the flux of unocculted light is greater underneath the flne-structure 
line (i.e, a broad C IV emission feature in the unocculted spectrum). 

Thus, the unocculted spectrum must consist of continuum, broad emission lines, and 
broad absorption lines similar to those of the nuclear spectrum itself. The most straightfor- 
ward explanation is that there is a scattering region that is outside of the region responsible 
for the D+E component, which reflects a fraction of the nuclear spectrum into our line 
of sight. This is a reasonable conclusion, given evidence for electron-scattered radiation in 
many Seyfert 2 galaxies (Moran ct al. 2000, and references therein) and the troughs of broad 
absorption-line (BAL) QSOs (Ogle et al. 1999). Additional evidence for electron scattering 
of the nuclear radiation in NGC 4151 comes from ground-based spectropolarimetry, which 
flnds signiflcant (~1%) polarization in the continuum and broad emission lines in relatively 
large (~3") apertures (Schmidt & Miller 1980; Axon et al. 1994; Martel 1998). Henceforth, 
we will refer to the unocculted emission in our echelle spectra as "scattered hght" . 

Since electron scattering is independent of wavelength, the spectrum of the scattered 
light in NGC 4151 should be the nuclear spectrum averaged over some time interval (except 
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for possible velocity shifts or broadening due to thermal or bulk motion of the scatterers). 
However, the extreme variability of the UV spectrum of NGC 4151 (Ulrich et al. 1991) and 
the rather loose limits on the location of the scattering region imposed by our aperture size 
(projected distances < 0".l or 6.4 pc from the central continuum source) make it difficult to 
determine what the time-averaged spectrum should look like. Given that NGC 4151 was in 
moderate to high states for at least 7 years prior to the STIS echelle observations (Weymann 
et al. 1997; Paper I) and evidence that the low-ionization absorption lines appear at only 
low continuum flux levels (Paper I), the most reasonable assumption is that the scattered 
spectrum resembles NGC 4151 in (at least) a moderate flux state and, hence, contains only 
the high-ionization broad absorption lines. 

We have adopted our STIS low-dispersion spectrum, which was obtained in 1998 at a 
high state and covers the entire UV region at high signal-to-noise (Nelson et al. 2000a; Paper 
I), as the basis of our scattered spectrum. In fitting the scattered hght profile, different scale 
factors are needed for different continuum regions and emission lines. For the continuum, 
we find that the best scale factors increase linearly from 0.12 to 0.20 over the 1150 - 3100 A 
region. We find scale factors of 0.12 for the high-ionization emission lines (e.g., C IV, N V, 
Si IV) and 0.4 for the low-ionization lines (e.g., Mg II, C III], the Fe II emission complex 
beginning at 2300 A). Nevertheless, there is observational justification for different scale 
factors. The time-averaged nuclear spectrum is hkely to reflect a lower state than our high- 
state spectrum, and the variability amplitudes of the low-ionization emission lines are much 
smaller than those of the continuum or high-ionization emission lines (Clavel et al. 1990). 
Also, the variability amphtude of the continuum in NGC 4151 (and other Seyfert 1 galaxies) 
decreases with wavelength in the UV (Crenshaw et al. 1996). Thus, larger scale factors are 
appropriate for the low-ionization emission lines and the long-wavelength continuum regions. 

We show our scattered-hght spectrum beneath the STIS echelle spectra in Figure 3 (a - 
c, the upper fit to the continuum and broad emission will be discussed later) . Qualitatively, 
the scattered-light spectrum provides a good match to the troughs of the broad absorption 
lines. In particular, the scattered spectrum no longer falls above the troughs of the high- 
ionization absorption lines, and it matches the troughs of absorption lines embedded in 
various emission lines reasonably well. For example, although the Mg II absorption appears 
not to be saturated at first glance, the scattered-hght spectrum shows that the A2796.3 hne 
is indeed saturated. Furthermore, tests show that inclusion of narrow Mg II emission in the 
nuclear spectrum, which should be present at the systemic redshift of NGC 4151 (Nelson et 
al. 2000a), would saturate the Mg II A2803.5 absorption as well. Although the scattered- 
light spectrum does not provide a perfect match to the troughs of the saturated lines, we 
conclude that that there is a significant flux of scattered light (~12% of the high-state FUV 
continuum) within the aperture, and that most of the broad absorption lines in the echelle 



-8- 



spectra must be highly saturated, indicating large ionic columns. In the observed low state, 
the percentage of light that is scattered is therefore very high (e.g. ~ 50% in the long 
wavelenth continuum). 

2.3. Absorption Measurements 

In this paper, we focus on the narrow kinematic components A, C, and E', and on 
the broad component D+E, since these components are likely to be photoionized (directly 
or indirectly) by the central continuum source. In addition, we present results on a new 
component mentioned in Paper I, which we will call D'; this component is shallow, broad 
(FWHM ^ 940 km s"^), highly blueshifted (-1680 km s~^), and only seen in the high 
ionization lines (N V, C IV, and Si IV). The D' component is not detected in previous STIS 
or GHRS spectra obtained at higher states but it is clearly seen as a flattening of the C IV 
absorption profile around 1544 A in the STIS echelle spectra (Paper I). Figure 4 shows the 
D' component in the N V A1238.8 absorption. This component was isolated by matching the 
D+E component with the profile derived from S II A1253.8 (Figure 2) and then subtracting 
this fit (in optical depth space) from the broad observed blend. 

We do not provide results on components B, F, and F' in this paper. Component B is 
from our Galaxy, and components F and F' are likely to arise from the interstellar medium 
and/or halo of NGC 4151 (Weymann et al. 1997); detailed analyses of these components 
will be provided in Danks et al. (2000). 

To determine the shape of the intrinsic continuum, we used a simple spline fit to regions 
that appeared free of absorption or emission. Similarly, several broad emission-line profiles 
were fit with cubic splines in regions unaffected by absorption, and an "average" broad profile 
was adopted (clearly this is a simplification, since the broad emission lines do not necessarily 
have the same profile in individual sources). This profile was then reproduced at the known 
positions of each broad line, and scaled in fiux until a suitable match was obtained. Our 
resulting estimate of the intrinsic continuum plus emission spectrum is given as the upper 
plot in Figure 3. 

Measurements of the narrow absorption components A, C, and E' are straightforward, 
since the absorption lines are relatively unsaturated and resolved (i.e., FWHM > 10 km s^^). 
Here we assume that these components cover the entire continuum plus broad emission 
regions. Since the residual intensities of the narrow components are large, a non-unity 
covering factor would have only a minor effect on the derived columns. Each absorption-line 
profile is converted to optical depth as a function of radial velocity, and integrated to obtain 
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the ionic column density (Crenshaw et al. 1999). Velocity centroids and widths (FWHM) 
are also determined from the optical depth profiles (see Table 1). 

Measurements of the broad absorption lines are more complicated. It is extremely 
difficult to separate Weymann et al.'s D and E components (contrary to our expectations 
in Paper I), because they are broad and blended together. The velocity centroids of all the 
isolated broad lines are very similar, which indicates that both components are present in 
each line, and it is therefore not possible to isolate the D or E component in a particular line. 
Since any differences in the observed profiles of different ions can be explained by saturation 
effects, contamination by other lines (particularly Fc II), or the scattered light spectrum, we 
can find no evidence that these components arc physically distinct. We will therefore assume 
that this is one component (D+E). We used a cubic spline to fit the D+E component in the 
S II A1253.8 line, as shown in Figure 2. 

Another complication is our finding that most of the broad (D+E) lines are heavily 
saturated. Although we have a quahtative idea of what the scattered spectrum looks like, 
we do not know the scattered flux beneath these lines exactly, and therefore it is not possible 
to determine their ionic column densities. The only relatively isolated lines that are not 
clearly saturated are S II A1253.8, Ni II A1317.2, and some of the metastable Fe II lines. 
We determined the velocity centroid and FWHM of D+E from these lines, and the column 
densities by integration of their optical depths. To obtain absolute lower limits on the D+E 
components of other isolated but saturated lines, we used optical depth profiles obtained 
by assuming no scattered light (i.c, a covering factor of 1.0). For doublets that arc blended 
(e.g., C IV), we used the S II A1253.8 optical depth profile to match the two lines assuming 
the proper ratio of oscillator strengths. 

We have retrieved the GHRS spectra of Weymann et al. (1997) from the HST archives 
to determine the column densities and velocities of the absorption components in a consistent 
fashion. We have also measured the absorption in a high-resolution STIS slitless spectrum 
of the C IV region observed on 1997 May 25 (Paper I). We list our measurements of the 
velocity centroid and FWHM of each component in the echelle spectra in Table 1. These 
values are consistent with those of Weymann et al. (1997), except that his component C is 
much broader (FWHM 140 km s~^) than ours. Our own measurements of the GHRS and 
STIS spectra show no evidence for changes in the velocities and widths of any component 
that are larger than those expected from measurement and wavelength calibration errors 
(typically ±10 km s~^). We present a comparison of the slitless (STISl), echelle (STIS2), 
and GHRS (epochs 1 -5) ionic column densities and lower limits in Table 2. 
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3. Modeling The Absorbers 

Photoionization models for this study were generated using the code CLOUDY90 (Fer- 
land et al. 1998). We have modeled the absorbers as matter-bounded slabs of atomic gas, 
irradiated by the ionizing continuum radiation emitted by the central source. As per con- 
vention, the models are parameterized in terms of the ionization parameter, U, the ratio of 
the density of photons with energies > 13.6 eV to the number density of hydrogen atoms 
(uij) at the illuminated face of the slab. Each separate kinematic component was modeled 
with its own set of initial conditions, i.e., U, n^, and the physical depth of the slab. As in 
earlier studies (e.g. Kriss 1998; Crenshaw & Kraemer 1999), a model is deemed successful 
when the predicted ionic column densities provide a good match (i.e., better than a factor 
of 2) to those observed. 

The SED for the ionizing radiation in NGC 4151 is difficult to determine, due to the 
combined effects of intrinsic absorption in the UV (Kriss et al. 1992, 1995) and X-ray (George 
et al. 1998, and references therein). Although the narrow emission line ratios are sensitive 
to the SED, Alexander et al. (1999) and Kraemer et al. (2000a) demonstrated that there is 
evidence for considerable absorption between the central source and the NLR. Hence, while 
both Alexander et al. and Kraemer et al. argued against the presence of a "Big Blue Bump" 
at energies > 13.6 eV, the narrow hues cannot be used to constrain the SED. As a result, we 
have assumed the simple "intrinsic" SED described in Kraemer et al. (2000a), specifically 
a series of power-laws of the form F,^ oc with a — 1 below 13.6 eV, a — 1.4 over the 
range 13.6 eV < hu < 1000 eV, and a = 0.5 above 1000 eV. 

As discussed in Paper I (and references therein), the UV continuum of NGC 4151 is 
highly variable. During the time of these observations, the extinction-corrected continuum 
flux at 1450 A was fuso = 1.1 x 10~^^ ergs cm~^ s~^ assuming an extinction E^.y = 
0.04 (Kriss et al. 1995). This is roughly 25% of the flux of NGC 4151 in a "high" state, e.g. 
during the GHRS epoch 5 observations (Weymann et al. 1997). In several cases, we have 
compared our models of the current ("low" state) absorbers with the high state absorbers 
observed with the GHRS; in order to do so, we increased the ionizing flux by a factor of 4, 
while, for simplicity, using the same SED, although there is some evidence for variations in 
the latter as a function of continuum flux (Edelson et al. 1996). 

Since La is saturated for most of the kinematic components, the ionization state and 
effective hydrogen column density (Ne// = N(H I) -|- N(H II)) of the absorbers is derived 
from the ionic columns of the heavy elements^° and, therefore, is sensitive to the elemental 



'We use N(XM) to denote ionic columns, where "X" is the atomic symbol and "M" is the ionization 
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abundances. Since there is no strong evidence for unusual abundances in the nucleus of 

NGC 4151 (although see Kraemcr ct al. 2000a), we have assumed roughly solar element 
abundances (cf. Grevesse & Anders 1989). They are, by number relative to H, as follows: 
He = 0.1, C = 3.4 X 10-^ N = 1.2 x 10-^ O = 6.8 x lO""^. Nc = 1.1 x 10~^, Mg = 3.3 x 
10-^ Al = 2.96 X 10-^ Si = 3.1 x 10-^ P = 3.73 x 10"^ S = 1.5 x IQ-^, Fe = 4.0 x 10-^ 
and Ni = 1.78 x 10~^. We have assumed that the gas is free of cosmic dust, hence all the 
elements are fully in the gas phase. 



4. Compcirison of the Observations and Model Results 

4.1. Component A 

GHRS spectra of Component A revealed absorption by C IV, Mg II, and Si II (Weymann 
et al. 1997). Kriss (1998) was able to match the Mg II and C IV columns with a single- 
component photoionization model, with U = 0.001 and Ne// = 2.0 x 10^^ cm^^. The STIS2 
data show absorption from the following ions: H I, Si II, Si III, Si IV, C II, C IV, and 
Mg II. We detected both C II A1334 and C II* A1336; the ratio of N(C II)/N(C IF) is ^ 
1.2, indicating that Ug (electron density) is ~ 100 cm~^ (Srianand & Petitjean 1999), which 
is consistent with the N(Si II)/N(Si IF) ratio, within the errors. After Shields & Hamann 
(1997), the C IV recombination timescale is t ~ (ug N(C IV)/N(C III) ac iv)~^- Using aciv 
(the radiative recombination rate) from ShuU Sz van Steenberg (1982), t is ~ 800 yrs. We 
measured a C IV column of (4.1 ± 0.3) x 10^^ cm~^, which is close to the earher GHRS 
results (see Table 2) as expected, given the long recombination time. The Mg II and Si II 
columns are identical to those from GHRS epoch 5, within the measurement errors. 

We were able to fit the ionic columns densities of A in the STIS2 spectra with a single 
component photoionization mode assuming U = 0.0012, and Ne// = 1.3 x 10^^ cm~^, similar 
to Kriss (1998). The model predictions are compared to the measured columns in Table 
3 (note that the model predicts total ionic columns, rather than separate predictions for 
ground and excited states of ions). The predictions matched the measured columns quite 
well, although N(Si III) is overpredicted, which may be due to the difficulty in determining 
the profile of the blue wing of La. 

As shown in Table 2, between GHRS epochs 4 and 5 there was a significant drop in the 
Mg II and Si II columns. Given the low densities derived for Component A, these changes 
are most likely the result of transverse motion, rather than a response to changes in the 



state. 
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ionizing continuum. We were able to match the GHRS epoch 4 results by adding a second 
component, with U = 0.0003 and Ne// = 4 x 10^^ cm~^; assuming the two components are 
at the same radial distance, the additional component would have a density of ~ 400 cm~^. 
The additional low ionization gas contributes no C IV absorption, as required by the lack 
of variation of N(C IV). This component must have a transverse velocity of ~ 2000 km s~^ 
in order to pass out of our line of sight to the broad line region (BLR) within 1.5 years, 
assuming the BLR of NGC 4151 is 4 light days in diameter (Clavel et al. 1990); this is not 
unreasonable, considering that the radial velocity of Component A is ~ 1588 km s~^. 

4.2. Component C 

Both C IV and Mg II absorption lines were identified with Component C in the CHRS 
spectra (Weymann et al. 1997). Based on these ionic columns, Kriss (1998) modeled the 
absorber as a single slab with U = 0.002 and N^ff — 1.25 x 10^^ cm~^. Our re-examination 
of the GHRS spectra revealed the presence of absorption by Si II A1526 in Epoch 1. In 
addition to these ions, the STIS2 spectra show absorption from Si III, Si IV, and C II. Based 
an the upper limit for C II* A1336, N(C II*)/N(C II) is - 0.2, which indicates Ue ~ 10 cm'^. 
At such low densities, Si II* A1533 should be weak or absent (Flannery et al. 1980), which is 
indeed the case. Also, the ionic columns should not respond to the changes in ionizing flux, 
which is true for Mg II, Si II, and C IV (Table 2). 

We modeled Component C as a single slab, with U= 0.0012, and Ng// = 1.0 x 10^^ 
cm^^. The model predictions and measured ionic columns for are listed in Table 4. Most of 
the ionic columns are well fit, albeit with overpredictions of Si III (see above) and Si II. 

4.3. Component D+E 

As noted in Section 2, we have not tried to deconvolve the Weymann et al.'s Components 
D and E. The cumulative absorption component (D+E) possesses strong high ionization lines, 
e.g. N V, C IV, and Si IV. In our STIS2 spectra, a number of lower ionization absorption lines 
are associated with Component D+E, the most striking of which are those from metastable 
levels of Fe II (Paper I), including those with lower levels up to 4.1 eV above the ground 
state. With the resolution of the STIS echelle, we have now established that the C III* A1175 
absorption, which arises from a level 6.5 eV above the ground state, is also associated with 
D+E, indicating densities of Ug > 10^ cm~^. 

The population of metastable states of C III and Fe II likely result from coUisional 
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excitation (although, for a discussion of the excitation of Fe II metastable states by UV 
radiation, see Verner et al. [1999]). The relative populations of the C III metastable, ^P", 
and ground, ^S^ states, as a function of Ue and temperature, have been calculated by Bromage 
et al. (1985) and Kriss et al. (1992). For a gas characterized by Ug = 10^'^ cm^^ and T = 
2 X 10"^K, a few percent of the C 111 will be in the excited state (while collision rates for 
temperatures below 2 x 10'^ K are unavailable at the time of this writing, it is probable that 
the relative population of the ^P" state will be smaller at lower temperatures). Under such 
conditions, at least the lowest 25 energy levels for Fe II (< 2.64 eV above the ground state) 
will be in statistical equihbrium (Verner et al. 1999). The ratio of Si IP/Si II cannot be 
used as a density indicator, since it is clear that these hues are saturated (see Section 2). 

As discussed in Section 2, the measurement of the ionic column densities for Component 
D+E is hampered by the complex profile of the underlying scattered component and the fact 
that the lines from the most abundant ions are saturated. The only isolated unsaturated 
lines are are S II A1254 (which has a weak oscillator strength, / = 0.01; Morton et al. 1988), 
Ni II A1317 (which has a low abundance; Section3), and some of the weaker Fe II metastable 
lines. We have, therefore, used the S II and Ni II columns to constrain the fraction of low 
ionization gas in the absorber. Although there are no weak high ionization lines in the STIS 
data, P V AA1117.98, 1128.01 has been detected in HUT Astro-1 and Astro-2 observations 
(Kriss et al. 1992; Kriss et al. 1995) and with the Berkeley sepctrometer aboard Orbiting and 
Retrievable Far and Extreme Ultraviolet Spectrometers (ORFEUS) - Shuttle Pallet Satellite 
(SPAS) II mission (Espey et al. 1998). Based on the equivalent widths determined by 
Espey et al., and the appearance that the hnes are unsaturated, N(P V) is 1.2 x 10^^ 
cm~^. Since the ionization potentials of P IV and C III arc similar (51.4 eV and 47.9 eV, 
respectively; Allen 1973), it is reasonable to expect that the ratio of ionic columns of P V 
and C IV will be approximately the same as that of the elemental abundances of phosphorus 
and carbon (i.e. 0.001; cf. Grevcssc & Anders 1989). Hence, at the time of ORFEUS-SPAS 
//observations, during which NGC 4151 was in a high flux state, N(C IV) must have been 
~ 10^^ cm~^. Although there are no far-UV observations contemporaneous with the STIS 
echelle observations, we will show that it is plausible that the columns of C IV and P V have 
remained relatively constant since the time of the ORFEUS-SPAS II mission. 

Based on these assumptions, we have modeled the Component D+E as a single slab, with 
U = 0.015, nn = 3.2 x 10*^ cm-^ and Ne// = 2.75 x 10^^ cm"^. The predicted and measured 
ionic columns are listed in Table 5. We ran two models, one assuming thermal broadening 
and the other a microturbulence velocity of 435 km s~^, and found no significant differences 
in the predicted ionic column densities (the results we discuss are from the former model). 
The model is successful in reproducing the Ni II and S II columns, within the measurement 
errors. In all other cases, with the exception of Al III, the model predictions exceed the 
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lower limits for the ionic column densities, which is a necessary, hut insufficient test of the 
model. The predicted average electron temperature is ~ 1.6 x lO^K, hence our assumption 
that fractional population of the state of C III is ~ 1% is reasonable. The predicted 
N(A1 III) is close to the derived upper limit, which may indicate that the Al III lines are 
not highly saturated, or that the aluminum may exceed solar, as has been suggested for 
BALQSOs (Shields 1997). 

As an additional test of the model, we ran a high state version, increasing U by a factor 
of 4 (see above), while holding Ngjj fixed. The results are listed in Table 5, alongside the 
low state model predictions. Notably, the Fc II and Si II columns arc small enough that 
their absorption lines would be weak or below the limit of detectability, which is precisely 
the case for the STISl (Paper I) and GHRS (Weymann et al. 1997) observations, which were 
all obtained at a high state. On the other hand, Mg II absorption would still be present; 
the predicted column density is similar to that seen in GHRS epoch 5 (which we determined 
to be 1.0 X 10^^ cm-2). The higher ionization lines, such as N V AA1238, 1241 and C IV 
AA 1548, 1551, would remain heavily saturated, and hence do not change appreciably in 
profile, as is the case (see Paper I). Most important, the C IV column density docs not 
change significantly, while N(P V) is close to the value derived from the ORFEUS - SPAS II 
observations. Therefore, given this model solution, it is possible to use the high state value 
of N(P V) to constrain Ng//. Although there may be other possible solutions, we will show 
in section 4.4 that the physical conditions in Component D' depend strongly on the column 
density and ionization state of Component D+E, and hence provide a further test. 



4.4. Component D' 

The appearance of Component D' was discussed in Section 2.3, and, as noted, it is 
present in N V, C IV, and Si IV (see Table 6). Since N(Si IV) is = 4.5 x 10^^ cm-^, there 
must be a significant column of C III, although no corresponding C III* A1175 is detected, 
indicating that the electron density is < 10^-^cm~^. Assuming that this component is directly 
ionized, we were unable to fit the measured ionic columns with a single-component model. 
Our best fit was a two component model, with the C IV and Si IV absorption arising in one 
component (U = 0.001; Ng// = 2 x 10^*^ cm^^), and the N V arising in more highly ionized 
gas (U > 0.1, Ng// > 1 X 10^° cm^^). However, if the ionizing fiux is increased by a factor 
of 4, the C IV column predicted by the lower ionization component increases, due to the 
ionization of the large C HI column; hence a directly ionized, multi-component model does 
not predict the absence of Component D' while NGC 4151 is in a high state. 

The coincidence of the appearance of D' at the same time as the low ionization lines 
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from D+E can be best explained if the former is screened by the latter. To demonstrate 
this, we generated a single-component model, with Ng// = 1 x 10^° cm^^, using the ionizing 
continuum filtered by D+E as the input spectrum (see Figure 5a). In order to reproduce the 
observed high ionization lines, assuming that these components are at the same approximate 
radial distance, the D' component must have U = 0.012, which yields a density of iih = 2.6 
X 10"^ cm~^. As shown in Table 6, the model predictions provide a good fit to the N V, C IV, 
and Si IV columns, and indicate that Si II and Mg II lines should be weak or absent (due to 
the small column densities and broadness of this feature). A second model was generated, 
using the filtered continuum from the high state model for Component D+E as the input 
(Figure 5b), while holding hh and Ng// fixed, which results in U = 4.15. The results are 
listed in the 3rd column of Table 6 and, as required by the observations, the predicted column 
densities are quite small. Therefore, the appearance of D while NGC 4151 was in its low 
state is the direct result of the higher EUV opacity of D+E . 



4.5. Component E' 

Component E', isolated kinematically by Weymann et al. (1997), is a low ionization 
absorber, as evidenced by the presence of Mg II A2803 and absence of C IV AA 1548,1551. 
Based on these constraints Kriss (1998) determined the ionization parameter for this com- 
ponent to be U < 0.0004, with ^eff ~ 10^^ - 10^*^ cm^^. In addition to Mg II, we have found 
absorption hues of O I, C II, Si II, and Fe II associated with this component, but, as with 
the GHRS results, no high ionization lines (Table 7). Interestingly, all of the ionic species 
present arise from states with ionization potentials below 13.6 eV. Among the more interest- 
ing features is the presence of absorption by the fine structure line O I* A1304 (in addition 
to the ground state line O I A1302). The ratio of N(0 I*)/N(0 I) is ^ 0.64, which is close 
to the ratio of the statistical weights of the two levels, indicating that they arc in collisional 
equilibrium. Based on the transition probabilities and collision strengths (Pradhan & Peng 
1995, and references therein). Component E' must have Ug > 10^ cm~^. This is consistent 
with the density indicated by the N(Si II*)/N(Si II) ratio, ~ 2, within the measurement 
errors, requiring Ug > 10^ (Flannery et al. 1980). 

We modeled E' as a single slab, with U = 0.0001, iih = 1 x 10® cm^"^, and Ng// = 1.5 x 
10^*^ cm^^. The model predictions are listed in Table 7, and show good agreement with the 
measured ionic columns. The predicted Fe II column is greater than that measured, which 
is evidence for a substantial population in metastable states, as one would expect at this 
density (Wampler et al. 1995) (the metastable Fe II lines were not detected in our spectra, 
due to their substantially lower oscillator strengths). A high state comparison model was 



-16- 



run, with U increased by a factor of 4. The predictions for N(Mg II) and N(Si II) were 6.7 
X 10^^ cm~^ and 1.3 x 10^^ cm~^, respectively, in good agreement with the GHRS epoch 4 
values (see Table 2). This indicates the increase in ionic column density of the low ionization 
lines between GHRS4 and STIS2 is likely due to a decrease in the ionizing flux. The columns 
for the high ionization lines remain small, e.g., N(C IV) = 7.4 x 10^^ cm~^, which would 
make them difficult to detect. 



5. Discussion 

5.1. The Geometry of the Circumnuclear Gas 

The presence of absorption from both ground and fine structure absorption lines of O I, 
Si II , and C II in several kinematic component, and the absorption from metastable C III and 
Fe II in Component D+E, allows us to constrain the densities of the absorbers. Combined 
with the evidence for screening of Component D' by Component D+E, and a significant 
contribution from an unoccultcd scattering region, we can begin to define the structure of 
the circumnuclear gas in NGC 4151. 

Based on the flux near the Lyman limit and our assumed SED, the number of ioniz- 
ing photons emitted by the central source is Q = 2 x 10^^ s~^. Prom our determination of 

their densities^^ and ionization parameters, wc place the absorbers at the following radial 
distances: Component C, 2.15 kpc; Component A, 681 pc; Component E', 23.6 pc, and 
Component D+E, 0.03 pc. As noted in the previous sections, we have assumed that Com- 
ponents A, C, and E' arc directly ionized by the central source, and hence not screened by 
intervening gas. However, it is entirely plausible that A and C arc screened by Component 
D+E. To test this, we generated models for these components using the high state filtered 
continuum (Figure 5), and found little change in the resulting column densities. The main 
effect is that these components would he at ~ 2/3 the radial distance determined from the 
direct ionization models. Due to its greater density. Component E' must respond to short 
timescale flux and opacity changes, and we were unable to model the small oberved changes 
in the ionic columns by assuming that this component is screened by D+E. Therefore, we 
suggest that E' is directly ionized, but not in our line-of-sight to the BLR and central source 
and, instead, is only occulting scattering emission^^ (which contributes 1/3 - 1/2 of the UV 



^^Note that the densities for Components D+E and E' arc lower hmits. 

^^If this is the case then the changes in the measured ionic columns over the last several years could at 
least partially be a result of an increased contribution of the scattered continuum to the total UV flux. 
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flux in the STIS2 data; see Section 2.2). Based on our model results, the absorbers decrease 
in density with increasing radial distance, faster than distance"^ up to Component A. The 
relative positions of the absorbers and the scattering region are shown in Figure 6. 

5.2. The Scatterer 

Based on our analysis of the scattered profile, we can begin to constrain the physical 
characteristics of the scatterer. Assuming isotropic scattering by free electrons (see Section 
2.2), for small electron scattering optical depths (r < 1), the reflected fraction of continuum 
radiation is ~ ^eFc^'thompson (the product of the column density of free electrons, the covering 
factor, and the Thompson cross-section, respectively). In order to scatter ~ 15% of the UV 
radiation (sec Section 2), the scatterer must have Ng ~ 2 x 10^'^ cm^^ (for a covering factor 
of unity). From STIS G430M spectra taken on May 15, 2000, we determined that the total 
broad + narrow H/3 flux from the central 0".l was ^ 5.8 x 10~^^ ergs s~^ cm~^ (Nelson et al. 
2000b), or Lho ~ 1.2 x 10^^ ergs s"^ (the distance to NGC 4151 is 13.3 Mpc, assuming H,, 
= 75 km s"^ Mpc~^). Since reverberation mapping indicates that the broad- line emission 
arises in relatively dense gas (cf. Clavel et al. 1990), which will not contribute large columns 
of free electrons, we may assume that the scatterer makes only a small contribution (i.e., 
< 10%) to the broad H/9 emission. Since the total H/? emission depends on the number 
of recombinations per unit volume, this provides an upper limit to the density. Assuming 
Component D+E is ~ 0.03 pc from the nucleus, and does not cover the scattering region, 
we may assume that the scatterer is coincident or lies further from the nucleus, e.g., 0.1 pc. 
Finally, this component cannot be opaque to the ionizing radiation, since there is strong 
evidence that the NLR is powered by the central source (Alexander et al. 1999; Nelson et 
al. 2000a; Kraemer et al. 2000a). With these constraints, the scatterer can be characterized 
with the following parameters'^: U > 4.5, Ng// = 2 x 10^^ cm~^, and uh < 1.0 x 10^ 
cm~^. Recent Chandra High Energy Transmission Grating spectra of NGC 4151 reveal an 
extended ('^ 1.6 Kpc), X-ray emission line region (Ogle et al. 2000), comprised of a two- 
phased medium: a hot (lO^K), collisionally ionized plasma, in which high ionization lines 
such as Fe XXV 1.79 A arise, and photoionized gas at temperatures < 4 x lO'^K, which give 
rise to emission lines and radiative recombination continua from lower ionization species. It 
is possible that scatterer is part of the hot medium, although Ogle et al. find little scattered 
continuum below 3 keV. 



^^The temperature at the illuminated face of this component is 2.2 x lO^K, hence it is not thermally stable 
(Krohk, McKee, & Tarter 1981) 



5.3. The X-ray absorber 



The scatterer model predicts large O VII, O VIII, Ne IX, and Ne X columns (9.7 x 
10^^ cm~^, 2.3 X 10^^ cm~^, 1.7 x 10^^ cm~^, and 3.5 x 10^^ cm~^, respectively), and, hence 
would appear as a "warm absorber" if viewed in our line-of-sight. This leads to the question 
whether the X-ray absorber and the scatterer are the same component, as suggested by 
Krolik & Kriss (1995) for Seyfert galaxies in general. This is particularly interesting since 
none of the UV components possesses sufficient X-ray opacity to account for the observed X- 
ray absorption. For example, Component D-|-E has the greatest column density and highest 
ionization parameter, but it transmits a considerable fraction of the intrinsic continuum at 
energies > 300 eV (see Figures 5a and 5b), which is clearly not observed in NGC 4151 
(George et al. 1998). Furthermore, as noted above, the scatterer must have a high covering 
factor, and hence, it is likely that it covers the central source. To address this, we examined 
the ASCA observations of NGC 4151 from 1993 November 5. Since NGC 4151 was in a 
higher flux state at this time than during the STIS2 observations, we can assume that the 
scatterer was in at least as high a state of ionization as the one we have modeled. However, 
the predicted ratios of N(0 VIII)/N(0 VII) and N(Ne X)/N(Ne IX) are too large to flt the 
aborption profile in the ASCA data, and the fit cannot be improved by varying the fraction of 
scattered light, column density of the scatterer, or the X-ray photon index (F). Instead, the 
best fit parameters are as follows: U = 0.86, Ng// = 4 x 10^^ cm^^, F = 1.5, and a scattered 
fraction of ^ 0.045 of the intrinsic flux (for the details of the analysis and fltting statistics, 
see George et al. 1998). The predicted ionic column densities for the X-ray absorber are 
hsted in Table 8. The model predicts large C IV and N V columns, and hence could explain 
some of the UV absorption. Nevertheless, N(P V) and N(Mg II) are too low, even for a high 
state period, hence we believe that this component in not responsible for the the bulk of the 
UV absorption. Furthermore, the C IV and N V absorption from the X-ray absorber could 
be buried in the deep D-l-E absorption trough, and hence undetectable, as suggested by Kriss 
et al. (1992). In summary, in addition to the UV absorbers and highly ionized scatterer, 
there is a third component of gas along our line-of-sight with an intermediate ionization 
parameter that produces the bulk of the X-ray absorption, similar to the multi-phased NLR 
suggested by Ogle et al. (2000). 

The decoupling of the X-ray and UV absorption presents a problem: if the absorbers 
occult the same continuum source, they must screen one another. We ran a scries of pho- 
toionization models to test the effects of screening. If the X-ray absorber is screened by 
the Component D-l-E, we found it will drop signiflcantly in ionization as the continuum flux 
drops and the screening gas becomes increasingly opaque to the ionizing radiation, similar 
to the case for Component D'. However, we flnd no evidence of the huge low ionization 
columns (e.g. N(Ni II) > 10^^ cm~^) which would arise from the X-ray absorber in such 
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a situation. Hence, in this scenario, the X-ray absorber must be of such low density that 
it does not respond on timescales of several years (i.e., uh < 10^ cm~^) which, in turn, 
requires that the short timescale variations in the X-ray absorption (Weaver et al. 1994; 
George et al. 1998) are due solely to transverse motion. On the other hand, if the X-ray 
absorber screens Component D-|-E from the ionizing source, as the ionizing flux dropped, 
we found that the UV absorber becomes far too neutral, unless the columns density of the 
X-ray absorber decreases by at least a factor of 4. 

Although we cannot rule out the possibility that the UV and X-ray absorbers screen one 
another, wc suggest another scenario. If the regions producing the UV and X-ray continua 
are not equally extended, each region could be covered by the respective absorber, i.e., the 
UV aborber covers the UV source and the BLR clouds, while the X-ray absorber covers only 
the X-ray source. Since the UV absorber must cover the BLR (~ 4 light days; Clavel et 
al. 1990), the simplest picture would have the X-ray absorber closer to the central source, 
occulting a smaller X-ray emitting region. The problem with this scenario is that the fraction 
of scattered light is apparently smaller in the X-ray than in the UV (~ 0.04 versus ~ 0.15 
at a high state). However, this could be the result of 1) a (relatively) anisotropic X-ray 
continuum, as compared to the UV or 2) an isotropic X-ray absorber inside the radius of 
the scattering region, which modifies the X-ray continuum before it reaches the scattering 
region. In the latter case, a larger scattered light contribution would be necessary to fit the 
ASCA data. Since our fit to the UV data requires the presence of absorption profiles in the 
scattered light, and there is evidence for the presence of an X-ray absorber inside the NLR 
of NGC 4151 (Kraemer et al. 2000a), it would not be surprising that the scattered X-ray 
continuum would show strong absorption features. 



6. Summary 

We have used medium resolution echelle spectra obtained with HST/STIS to study the 
physical conditions in the intrinsic UV absorbers in the Seyfert 1 galaxy, NGC 4151, while 
it was in a recent low continuum state. The present paper is a detailed study following 
our observational results from Paper I. We have determined the following regarding the 
circumnuclear gas in NGC 4151. 

1. Each of the kinematic components can be modeled as a single slab of atomic gas, 
of roughly solar elemental abundances, ionized by the EUV - X-ray continuum radiation 
emitted by the central source. The combined kinematic component D+E, which we have not 
attempted to deconvolve, has the largest column density (2.75 x 10^^ cm^^) and ionization 
parameter (U = 0.015) of the absorbers. The C HI and Fe II metastable lines arise in this 
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component, indicating it is also the densest (n^^ > 3 x 10^ cm~^) and closest to the ionizing 
source. Component D', which was not detected in earlier high state spectra, appears to 
be screened from the nuclear source by Component D+E; it becomes too highly ionized to 
detect as the opacity of D+E decreases when NGC 4151 enters a high flux state. Our model 
predictions match the observed changes in the ionic columns; specifically that the denser 
components D+E, D', and E' responded to the decrease in continuum flux since the time of 
the GHRS spectra, while the low density components A and C, have not (as a result of their 
long recombination timescales). None of the UV absorbers has sufficient X-ray opacity to 
produce the observed X-ray absorption. 

2. We have been able to constrain the densities of the UV absorbers based on the ratios 
of fine structure to ground state C II, O I, and Si II lines, and the presence of fines from 
metastable states of C III and Fe II. Prom these densities, and the ionization parameters 
determined from the fits to the ionic column densities, we can determine the relative radial 
positions of the individual absorbers, as shown in Figure 6. Based on these results, it is 
apparent that the UV absorbers extend over a large range in radial distances, and decrease 
in density with distance. 

Weymann et al. (1997) found that the radial velocities of Components A and C had 
changed little between GHRS epochs 2 and 5. From this, they determined the limit for 
steady radial acceleration for Component A to be ~ 1 x 10^"^ cm s~^. Assuming that 
the cloud is optically thin, and a mean observed continuum flux for NGC 4151, this limit 
requires that Component A lies at a radial distance of > 0.7 pc. Since the distances we have 
determined for Components A, C, and E' are all significantly greater than this limit, we find 
no apparent confiicts with the Weymann et al.'s conclusions. Since Component D+E fies 
much closer to the nuclear source, one might expect that it should experience a noticeable 
acceleration. However, since D+E is optically thick in the Lyman continuum and has a large 
column density, it is probable that radiative acceleration will be less efficient (Williams 1972; 
Mathews 1974). Additionally, the evidence for transverse motion in the UV and X-ray is 
evidence for more complex dynamics, which we intend to address in a subsequent paper. 

3. Based on the residual fiux in the troughs of staturated absorption lines associated 
with Component D+E, we estimate that at least 15% of the high-state continuum emission, 
broad emission lines, and absorption lines are scattered into our line of sight, from a plasma 
that extends outside the solid angle subtended by the absorber (in Figure 6 we show a 
possible location for the scatterer). For a covering factor of unity, the scatterer must have 
a free electron column density of Ng// ~ 2 x 10^^ cm~^. Based on the constraints on 
the total H/? fiux, the covering factor, and the EUV opacity of the scatterer, we argue 
that it must be highly ionized (U > 4.5) and, hence, cannot produce the observed X-ray 
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absorption. Instead, a smaller column of lower ionization gas (Ng// = 4 x 10^^ cm~^, U = 
0.86) must also be present. There are several possibilities regarding the relative position of 
this additional component with respect to the UV absorbers and scatterer. Also, the reason 
for the apparent differences in the fractions of UV and X-ray scattered light remains an open 
question, which can be better addressed through simultaneous UV (HST/STIS) and X-ray 
{Chandra) observations of NGC 4151. 

S.B.K. and D.M.C. acknowledge support from NASA grant NAG5-4103. We thank Fred 
Bruhweiler, Tahir Yaqoob, and Jane Turner for enlightening discussions. 
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Fig. 1. — STIS echelle spectra of the nucleus of NGC 4151; the spectra have been smoothed 
by a boxcar function of 7 pixels for display purposes. Identifications of the broad absorption 
lines associated with D+E, except for Fe II, are given above the spectra. Fe II multiplets 
are plotted below the spectra, with the length of each vertical line representing the gf value 
from Silvas & Bruhweiler (2000). 

Fig. 2. — Plots of the Si IV A1393.8 and S II A1253.8 lines as functions of radial velocity 
(with respect to the nucleus of NGC 4151 at a rcdshift of 0.00332). The spectra have been 
normalized by dividing by the continuum plus emission fit (Figure 3). The expected positions 
of the absorption components arc plotted as vertical dotted lines. The smooth line in the 
lower plot is our fit to the D+E component of S II; note there may be a slight contamination 
at —100 to +300 km s~^ due to an Fe II line from multiplet UV CW. Component B is 
Galactic, and Components F and F' are most likely from the interstellar medium or halo of 
NGC 4151. 

Fig. 3. — STIS echelle spectra of NGC 4151. The lower plots give our estimate of the 
scattered light spectrum (see text). The upper plots give our continuum plus broad-emission 
fits to the echelle spectra. 

Fig. 4. — Plot of the N V A1238.8 absorption as a function of radial velocity. The upper 
smooth plot gives our fit to the D' component, the lower smooth plot gives the sum of the 
D' and D+E components. The location of N V emission from the narrow-line region (NLR) 
is noted. The scattered component has not been subtracted in this plot. 

Fig. 5. — The incident (sohd fine) and transmitted (dotted fine) continua for the models of 

Component D+E: a) the low flux state model and b) the high flux state model (see text. 
Section 4.4). The flux values for the incident continuum are in log(z/F,y) at the face of the 
slab (in units of ergs s~^ cm~^). 

Fig. 6. — Schematic diagram of the nuclear region of NGC 4151, showing the relative po- 
sitions of the central source, the emission-line bicone, and the UV absorbers. The viewing 
angle relative to the bicone is given by the kinematic model by Crenshaw et al. (2000a). The 
dotted line represents a possible radial position and distribution of the scattering medium, 
although the only real constraint is that it is not completely covered by D+E in projection. 
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Table 1. Absorption Components in NGC 4151 



Component 


Centroid 


FWHM 




(km s-i) 


(km s-i) 


A 


-1588 (±8) 


36 (±4) 


B 


-1025 (±5) 


a 


C 


-858 (±8) 


27 (±3) 


D+E 


-491 (±8) 


435 (±12) 


D' 


--1680 (±85) 


940 (±60) 


E' 


-215 (±11) 


59 (±8) 


F 


-16 (±4) 


34 (±4) 


F' 


+34 (±4) 


15 (±1) 



'^sensitive to saturation level, see Danks 
et al. (2000). 
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Table 2. Comparison of Component Ionic Columns (x 10^^ cm ^) from STIS and GHRS 

Spectra" 







VjrXirLoZ 






^xlxVoO 




o J- loz 




22 06 92'' 


04 07 92 


03 01 94 


28 10 94 


11 03 96 

J- J- t\JXJ »>J\J 


25 05 97 


15 07 99 


Line 






Comp. A 










Si II A1526 


1.22 










<0.20 


0.46"= 




(±0.22) 












(±0.08) 


Si II* A1533 


1.17 


0.97 


1.05 


1.27 


0.51 


<0.20 


0.53*^ 




(±0.22) 


(±0.24) 


(±0.19) 


(±0.28) 


(±0.16) 




(±0.08) 


C IV A1548 




2.83 


3.62 


3.90 


3.84 


3.13 


4.11 






(±0.31) 


(±0.34) 


(±0.17) 


(±0.22) 


(±0.18) 


(±0.31) 


Mg II A2796 






0.46 


0.48 


0.29 




0.27 








(±0.03) 


(±0.04) 


(±0.03) 




(±0.04) 


Line 






Comp. C 










Si II A1526 


0.29 










0.32 


0.46^^ 




(±0.13) 










(±0.12) 


(±0.06) 


Si II* A1533 


<0.20 


<0.20 


<0.20 


<0.20 


<0.20 


<0.20 


<0.20 


C IV A1548 




2.21 


2.65 


2.79 


2.67 


2.50 


2.95 






(±0.37) 


(±0.32) 


(±0.41) 


(±0.25) 


(±0.43) 


(±0.40) 


Mg II A2796 






0.09 


0.11 


0.12 




0.10 








(±0.03) 


(±0.02) 


(±0.03) 




(±0.03) 


Line 






Comp. E' 










Si II A1526 


<0.20 


<0.20 


<0.20 


0.34 


<0.20 


<0.20 


1.37 










(±0.11) 






(±0.19) 


Si 11* A1533 


<0.20 


<0.20 


<0.20 


0.63 


<0.20 


<0.20 


2.45 










(±0.17) 






(±0.31) 


C IV A1548 


<0.30 


<0.30 


<0.30 


<0.30 


<0.30 


<0.30 


<0.30 


Mg II A2803 






0.61 


0.68 


0.64 




1.90 








(±0.09) 


(±0.06) 


(±0.06) 




(±0.13) 



''blank entry: not in wavelength coverage. Values in parentheses are the measurement errors 
for the previous line 

^Observation date (UT) in Day.Month.Year. 
'^from Si II A1260.4. 
<ifrom Si 11* A1264.7. 
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Table 3. Predicted and Observed Ionic Columns (x 10^^ cm ^) for Component A 



Ion (Line) Model" Prediction Measured Column 



H I (L-a) 


8.12E+02 


> l.OE+Ol'' 


C II (A1334) 




3.43 (±0.53) 


C II* (A1336) 




4.19 (±0.46) 


C II (total) 


7.42 


7.62 (±0.70) 


C III 


3.25E+01 




C IV (A1548) 


3.84 


4.11 (±0.31) 


Mg II (A2796) 


0.21 


0.27 (±0.03) 


Mg II (A2803) 


0.21 


0.29 (±0.05) 


Si II (A1260) 




0.46 (±0.08) 


Si II* (A1265) 




0.53 (±0.15) 


Si II (total) 


0.93 


0.99 (±0.17) 


Si III (A1206) 


2.09 


0.55 (±0.07) 


Si IV (A1393) 


0.80 


0.98 (±0.12) 


Si IV (A1403) 


0.80 


1.12 (±0.22) 


Fe II 


0.02 





"^U = 0.0012, Ne// = 1.3E+18 cm-^. 
''Saturated. 
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Table 4. Predicted and Observed Ionic Columns (x 10^^ cm ^) for Component C 



Ion (Line) 


ModeP Prediction 


Measured Column 


H I (L-a) 


6.40E+03 




C II (A1334) 


5.84 


4.71 (±0.36) 


C III 


2.50E+01 




C IV (A1548) 


2.88 


2.95 (±0.40) 


C IV (A1551) 


2.88 


3.20 (±0.44) 


Mg II (A2796) 


0.17 


0.07 (±0.01) 


Mg II (A2803) 


0.17 


0.10 (±0.03) 


Si II (A1260) 


0.76 


0.30 (±0.06) 


Si III (A1206) 


1.57 


0.58 (±0.05) 


Si IV (A1393) 


0.61 


0.51 (±0.07) 


Si IV (A1403) 


0.61 


0.70 (±0.17) 



= 0.0012, Ne// = l.OE+18 cm-2. 
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Table 5. Predicted and Observed Ionic Columns (x 10^^ cm ^) for Component D + E 



Ion (Line") Model'' Prediction Highstate'^ Prediction Measured Column 



H I (L-a) 


6.4E+06 


1.6E+04 




C II 


2.0E+04 


1.43E+01 




C III* A1175 


3.9E+02'' 


5.0E+01 


> 3.8E+01 


C IV 


2.5E+04 


3.1E+04 


> 4.3E+02 


N V 


1.5E+03 


7.9E+03 


> 2.7E+02 


I 


2.9E+04 






Mg II 


1.6E+03 


6.3 


> 3.4E+01 


Al III 


l.lE+02 


3.4 


> 1.2E+02 


Si II (A1526) 






> 1.3E+02 


Si II* (A1533) 






> 8.8E+01 


Si II (total) 


2.1E+03 


0.21 


> 2.2E+02 


Si III (A1206) 


1.7E+03 


2.7E+01 


> 1.6E+02 


Si IV 


1.9E+03 


2.1E+02 


> l.lE+02 


P V 


7.7 


6.5 


1.2E+01^ 


S II 


6.3E+02 




l.OE+03 (±5, 


Fe II (A1608) 


2.3E+03 




> 1.7E+02^ 


Ni II 


l.lE+02 




8.5E+01 (±2, 



.OE+02) 
.2E+01) 



'^when no wavelength is listed, the measured column is derived from a simultaneous 
fit to two or more lines (see text). 

= 0.015, Ne// = 2.75E+21 cm'^. 
= 0.060, Ne// = 2.75E+21 cm'^. 

'^assuming NmetastaUe/'^ ground = 0.01 (sCC text). 

^value from ORFEUS-SPAS 11 ohseiyation (Espey at al. 1998). 
'determined from resonance lines. 



Table 6. Predicted and Observed Ionic Columns (x 10^^ cm ^) for Component D' 



Ion (Line") 


Model* Prediction 


Highstate'^ Prediction 


Measured Column 


H I (L-a) 


3.34E+04 


3.33 




C II 


8.54E+01 






C III 


3.83E+02 






C IV 


2.58E+02 


0.19 


1.20E+02 (±2.40E+01) 


N V 


1.70E+02 


0.97 


2.20E+02 (±4.40E+01) 


Mg II 


1.30 






Si II 


3.26 






Si III 


2.15E+01 






Si IV 


2.82E+01 




4.50E+01 (±9.00) 



'^when no wavelength is listed, the measured column is derived from a simultaneous 

fit to two or more lines (see text). 

'^screened by Component D+E (see text), U = 0.012, Ne// = l.OE+20 cm"2_ 

'^screened by Highstate Component D+E (see text), U = 4.15, Ne// = l.OE+20 
cm~^. 
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Table 7. Predicted and Observed Ionic Columns (x 10^^ cm ^) for Component E' 



Ion (Line) 


Model Prediction" 


Measured Column 


H I (L-a) 


8.83E+03 




C II (A1334) 




9.73 (±3.18) 


C ir (A1335) 




1.23E+01 (±3.11) 


C 11 (total) 


3.56E+01 


2.20E+01 (±4.45) 


C III 


1.46E+01 




I (A1302) 




5.47 (±2.54) 


I* (A1304) 




3.48 (±1.54) 


I (total) 


6.26 


8.95 (±2.97) 


Mg II (A2803) 


2.03 


1.90 (±0.13) 


Si II (A1526) 




1.37 (±0.19) 


Si ir (A1533) 




2.45 (±0.31) 


Si II (total) 


3.75 


3.82 (±0.36) 


Si III (A1206) 


0.86 




Fe II (A2383) 


2.31 


0.50 (±0.11) 


Fe II (A2600) 


2.31 


0.31 (±0.13) 



•^U = 0.0001, Ne// = 1.5E+18 cm-2. 




Fig. la. 




Fig. lb. 




Fig. Ic. 



Flux (10"'' ergs s"' cm"' A"') 




Fig. Id. 
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Fig. If. 
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Fie. 2. 



-40- 




-41 - 




-42 - 



Flux (10"^' ergs s"' cm"' A"^) 
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Fig. 3c. 




Fig. 4. 
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Log z/F 
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Fig. 5a. 
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Fig. 5b. 




Fig. 6. 



TABLE 8 

Predicted Ionic Column Densities (cm~^) from X-ray Absorber Model" 



Element 


I 

Y 

A- 


II 

YT 




III 

YTT 
Aii 




IV 

YTTT 




V 

Ai V 


VI 

A V 




VII 

A V J 




VIII 
XVII 


IX 

YA/TTT 
A V 111 


H 


9.5E+16 


4.0E-f 


22 


























He 


3.2E+14 


6.6E+ 


19 


3.9E4 


21 






















C 




1.0E+ 


12 


5.1E4 


15 


4.0E4 


17 


6.1E+18 


5.3E4 


18 


1.8E4 


18 








N 




1.5E+ 


11 


1.5E4 


15 


1.1E4 


17 


1.9E+17 


3.2E4 


18 


1.2E4 


18 


l.oE4 


17 









4.9E+ 


12 


7.2E4 


16 


5.6E4 


17 


7.5E+17 


3.5E4 


18 


1.9E4 


19 


2.7E4 


18 


2.0E+17 


Ne 




4.1E-f 


12 


3.6E4 


16 


3.2E4 


16 


1.8E+17 


1.1E4 


18 


8.5E4 


17 


8.9E4 


17 


1.3E+18 




9.9E+16 


2.7E-f 


15 


























Mg 




6.2E+11 


2.8E4 


15 


2.9E4 


15 


2.8E+15 


1.8E4 


16 


9.7E4 


16 


3.1E4 


17 


4.9E+17 




2.8E+17 


1.1E+ 


17 


3.9E4 


15 


4.7E+13 


















Al 




2.6E-f 


11 


1.7E4 


14 


1.7E4 


14 


1.1E4-15 


7.1E4 


15 


2.3E4 


16 


3.8E4 


16 


3.3E+16 




1.2E+16 


3.3E+15 


8.6E4 


13 


6.9E4 


11 


















Si 








1 7E-I- 


11 


2.1E4 


13 




l.lE+16 


fi 4F,-I- 


1 fi 

J-U 


2.3E4 


17 






3.3E+17 


1.8E+ 


17 


4.1E4 


16 


9.6E4 


15 


1.7E+14 


8.7E4 


11 












TO 
C 








1.3E4 


10 


2.7E4 


12 


1.9E+13 


1.1E4 


14 


0.UE4 


1 A 
14 


2.6E4 


15 


4.4E+15 




3.6E-I-15 


2.3E+ 


15 


8.3E4 


14 


1.9E+14 


2.9E+13 


3.6E+11 












s 








4.3E4 


11 


7.0E4 


13 


8.4E-I-14 


4.0E4 


15 


2.6E4 


16 


9.1E4 


16 


1.6E+17 




1.5E+17 


1.1E+ 


17 


4.9E4 


16 


1.2E4 


16 


1.6E+15 


2.4E4 


14 


2.1E4 


12 








Fe 








1.5E4 


10 


3.0E+14 


1.3E+15 


3.8E4 


15 


9.6E4 


15 


3.9E4 


16 


1.6E+17 




2.4E+17 


2.8E+ 


17 


3.2E4 


17 


2.5E4 


17 


1.6E+17 


9.3E4 


16 


2.6E4 


16 


1.4E4 


16 


1.3E+15 


Ni 












1.8E4 


12 


2.1E+12 


2.8E4 


13 


1.6E4 


14 


8.8E+14 


4.9E+15 




l.lE+16 


1.86E 


+ 16 


1.4E4 


16 


1.1E4 


16 


7.1E+15 


3.8E4 


15 


1.6E4 


15 


5.3E4 


14 


1.5E+14 


'^U = 0.8 


6, N,ff - 4 


X 1022 


cm" 


2 

























